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Abstract. The reaction pp^dK + K° has been investigated at excess energies Q= 47A and 104.7 MeV 
above the K + K° threshold at COSY- Jiilich. Coincident dK + pairs were detected with the ANKE spec- 
^ \ trometer, and subsequently ~ 2000 events with a missing K° invariant-mass were identified, which fully 

populate the Dalitz plot. The joint analysis of invariant-mass and angular distributions reveals s-wave 
1 1 ' dominance between the two kaons, in conjunction with a p-wave between the deuteron and the kaon pair, 

i.e. KK production via the a f (980) channel. Integration of the differential distributions yields total cross 
sections of a(pp—^dK + K )=(38±2 s t a t±U sys t) nb and (190±4 a t a t±39 ay st) nb for the low and high Q value, 
respectively. 
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PACS. 25.10.+S Nuclear reactions involving few-nucleon systems - 13.75.-n Hadron-induced low- and 
intermediate-energy reactions and scattering (energy < 10 GeV) 

1 Introduction Experimentally more states with quantum numbers 

J P =Q + have been identified than would fit into a single 
SU(3) scalar nonet: the / (600) (or er), / (980), /o(1370), 
QCD is the fundamental theory of Strong Interactions. / (1500) and / (1710) with 7=0, the putative k(800) and 
How quarks and gluons are bound into hadrons is as yet an the 7f*(1430) (1=1/2), as well as the a (980) and a (1450) 
unsolved strong-coupling problem. Though QCD can be (1=1) [3]. Consequently, the a // (980) have been asso- 
treated explicitly in this regime using lattice techniques [1], dated with crypto-exotic states like KK molecules [4] or 
these are not in a state to make quantitative statements compact qq-qq states [5] . It has even been suggested that 
about the light scalar mesons. Alternatively, QCD-inspired a complete nonet of four-quark states might exist with 
models, which employ effective degrees of freedom, can be masses below 1.0 GeV/c 2 [61. 
used. The constituent quark model is one of the most suc- 
cessful in this respect (see e.g. Ref. [2]). This approach The first clear observation of the isovector a (980) res- 
inherently treats the lightest scalar resonances a //o(980) onance was achieved in K~p interactions [7]. It has also 
as conventional qq states. been seen in pp annihilations [8], in ir~p collisions [9], 

and in 77 interactions [10]. Experiments on radiative <jy 

a Present address: Physikalisches Institut, Universitat Bonn, decays [11,12] have been analysed in terms of a //o pro- 

Nussallee 12, 53115 Bonn, Germany duction in the decay chain — * 7ao//o — * 77r°ry/7r 7r°. In 

b m.buescher@fz-juelich.de pp collisions the ao(980) resonance has been measured at 

c deceased p v — 450 GeV/c via /i(1285) — ► ag7r T decays [13] and 

d Present address: Research Center for Nuclear Physics, Os- in inclusive measurements of the pp — > dX + reaction at 

aka University, Ibaraki, Osaka 567-0047, Japan p p = 3.8, 4.5, and 6.3 GeV/c [14]. 
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Despite these many experimental investigations, basic 
properties and even the nature of the a (980) resonance 
are still far from being established (see e.g. Refs. [15,16]). 
The Particle Data Group quotes a mass of m ao — (984. 7 ± 
1.2) MeV/c 2 and a width of r a „ = (50 - 100) MeV/c 2 [3]. 
The main decay channels, ttt] and KK, are denoted as 
"dominant" and "seen", respectively. The corresponding 
coupling constants g wv and g K x differ significantly for the 
different data sets and analyses [16]. 

Therefore, an experimental programme has been started 
at the Cooler Synchrotron COSY-Julich [17] aimed at ex- 
clusive data on the ao//o(980) production from pp, pn, pd 
and dd interactions at energies close to the KK thresh- 
old [18]. The final goal of these investigations is the extrac- 
tion of the ao//o-mixing amplitude, a quantity which is 
believed to shed light on the nature of these resonances [19, 
20]. As a first step the reactions pp —* dK + K° [21] and 
pp — > d-K + r] [22] have been measured in parallel at the 
ANKE spectrometer [23] for T p =2.65 GeV, correspond- 
ing to an_ excess energy of Q=47.4 MeV with respect to 
the K + K° threshold. The data for the strangeness de- 
cay channel — which are almost background free — in- 
dicate that more than 80% of the kaon pairs are pro- 
duced in a relative s-wave, corresponding to the chan- 
nel [21]. On the other hand, the ir + r/ signal sits on top 
of a strong but smooth background of multi-pion produc- 
tion which, together with the small acceptance of ANKE 
for this channel, makes the interpretation of the sig- 
nal model-dependent [22] . However, the obtained branch- 
ing ratio cr(pp^d(K + K ) s ^ wavc )/<j(pp^da+^diT + ri) = 
0.029 ± 0.008 sta t ± 0.009 sys [22] is in line with values from 
literature [3]. 

In this paper we report on a refined analysis of the 
pp — > dK + K° data at T p =2.65 GeV as well as on new re- 
sults from a second measurement at higher beam energy 
(T p = 2.83 GeV, corresponding to Q = 104.7 MeV). The 
procedures for event identification and acceptance correc- 
tion at the lower energy have been described in our pre- 
vious publication [21]. 



2 Measurement of pp — > dK+K° events with 
ANKE 

2.1 Experimental setup 

ANKE is a magnetic spectrometer located in one of the 
straight sections of COSY and comprises three dipole mag- 
nets, Dl - D3 [23]. Dl deflects the circulating COSY beam 
onto the target in front of D2, and D3 bends it back into 
the nominal orbit. The C-shaped spectrometer dipole D2 
separates forward-going reaction products from the COSY 
beam and allows one to determine their emission angles 
and momenta. The angular acceptance of ANKE covers 
|$h| < 10° horizontally and \{) v \ < 3° vertically for the 
detected deuterons (p d > 1300 MeV/c), and < 12° 
and < 3.5° for the K + mesons. 



A cluster-jet target [24] of hydrogen molecules, placed 
between Dl and D2, has been used, providing areal den- 
sities of up to ^5xl0 14 cm~ 2 . The luminosity has been 
measured with high statistical accuracy using pp elastic 
scattering, recorded simultaneously with the dK + data. 
Protons with = 5.5° — 9° have been selected, since the 
ANKE acceptance changes smoothly in this angular range 
and the elastic peak is easily distinguished from back- 
ground in the momentum distribution. The average lu- 
minosity during the measurements with up to ~ 4 x 10 10 
stored protons in the COSY ring has been determined as 
L = (1.7±0.4 syst )xl0 31 s _1 cm -2 , corresponding to an in- 
tegrated value of L; nt = 7.5 pb -1 . 

2.2 Event selection for Q = 104.7 MeV 

Two charged particles, K + and d, have been detected in 
coincidence. Positively charged kaons are identified in the 
side detection system (SD) [23,25] of ANKE by a time- 
of- flight (TOF) measurement. The TOF-start counters, 
consisting of one layer of 23 scintillation counters, have 
been mounted next to the large exit window of the vac- 
uum chamber in D2. Kaons from decay with momenta 
Pk+ = 390 — 625 MeV/c have been stopped in range 
telescopes, located along the focal surface of D2. These 
telescopes comprise TOF-stop counters and provide ad- 
ditional kaon- us-background discrimination by means of 
energy-loss (AE) measurements [25]. At T p — 2.83 GeV, 
kaons with p K + — (625 — 1000) MeV/c have been de- 
tected in a different part of the SD, consisting of one layer 
of 6 scintillation counters for TOF-stop ( "sidewall coun- 
ters"). Two multi-wire proportional chambers (MWPCs) 
positioned between the TOF-start and -stop counters al- 
low one to deduce the ejectile momenta and to suppress 
background from secondary scattering [21,22]. 

Fast particles produced in coincidence with the K + 
mesons as well as elastically scattered protons have been 
detected in the ANKE forward-detection system (FD) [26] 
containing two layers of scintillation counters for TOF and 
AE measurements. In addition there are three MWPCs, 
each with two sensitive planes, exploited for momentum 
reconstruction and background suppression [21,23]. Two 
bands of protons and deuterons are distinguished in the 
time difference between the detection of a i4T + -meson in 
one of the TOF-stop counters of the SD and a particle in 
the FD as a function of the FD particle momentum, see 
Fig. la. The deuterons are selected with the cut indicated 
by the lines, plus the energy-loss information from the FD 
scintillation counters. In Fig. lb the missing-mass distri- 
bution m{pp 1 dK + ) for the selected pp^dK + X_ events is 
presented. The missing particle X must be a K°, due to 
charge and strangeness conservation. The measured dK + 
missing-mass distribution peaks around m—m^-o, reflect- 
ing the clean particle identification at ANKE. 

About 2300 events are accepted as dK + K° candidates 
for further analysis (unshaded peak area of the histogram 
in Fig. lb)). The remaining background from misidentified 
particles is (13±2) %. The shape of this background in the 
differential spectra discussed below has been determined 
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Fig. 1. a) Time difference between the fast forward-going par- 
ticles in layer 1 of the FD scintillators and the K + -mesons 
vs. the momentum of the forward particle. The lines indi- 
cate the selection for deuteron identification, b) Missing-mass 
m(pp, dK + ) distribution of the pp—>dK + X events. The shaded 
areas indicate the events used for background subtraction, the 
solid line shows the background distribution under the K° peak 
obtained from a polynomial fit. 



value of mgo = 497.6 MeV/c 2 on an event-by-event ba- 
sis, varying the momentum components of the detected 
K + and d within their resolutions. As a result of the fit, 
the deuteron missing-mass (i.e. the invariant K + K° mass) 
resolution improves from 8m K+R0 = (35 — 3) MeV/c 2 over 
the range (991 - 1096) MeV/c 2 to Sm K+Ro < 10MeV/c 2 
in the full range with minimum values of ~ 3MeV/c 2 at 
the kinematic limits. Due to the fact that the p z resolu- 
tion (z being the beam direction) for deuterons is approx- 
imately a factor five worse than for all other variables, 
the fit procedure does not significantly improve the K + 
missing-mass and angular resolutions: Sm dR o ~ 5 MeV/c 2 
in the full range, <5[cos(#)] ~ 0.2 for all angular spectra. 

The same fit procedure has also been applied to the 
previously published data at T p =2.65 GeV [21], and im- 
proves 5m K +Ko from (8 - 1) MeV/c 2 over the range (991 
- 1038) MeV/c 2 to 8m K + R0 < 3MeV/c 2 in the full mass 
range. 



2.4 Non-acceptance-corrected Dalitz plot 

Figure 2 shows the distribution of the kinematically fit- 
ted dK + K° events in the Dalitz plot for both Q val- 
ues. It is observed that the kinematically allowed region 
is fully covered by the ANKE acceptance. For compar- 
ison the simulated population of the Dalitz plot is also 
shown for the case of phase-space-distributed events. The 
total ANKE acceptance for these dK + K° events is 2.1% 
at Q = 47 .4 MeV and 0.8% at 104.7 MeV. Due to the 
limited number of counts we present in the following only 
one-dimensional distributions. These also have the advan- 
tage of carrying additional information about the transi- 
tion matrix, as shown in Sect. 2.5. 



and subsequently subtracted by selecting events outside 
the K° peak (shaded areas in Fig. lb). 

The K + tracking efficiency in the side MWPCs and 
the efficiency of the AE cut have been determined using 
simultaneously recorded pp^pK + A events, which, due to 
the significantly larger cross section, can be identified by 
TOF criteria and a momentum cut for protons in the FD 
only. The efficiency of the track reconstruction varies from 
96% for the telescopes to 76% for the sidewall counters. 
The efficiency of the AE cut has been determined for each 
telescope, with an average value of 53%, and for the side- 
wall counters, ranging from 87% to 74%. The efficiency of 
the FD AE criterion for deuterons has been deduced from 
the number of K° events in the peak of Fig. lb before and 
after this cut. The efficiency of the FD scintillators and all 
TOF criteria is larger than 99%. The data have been cor- 
rected for all efficiencies on an event- by- event basis. 



2.3 Kinematic fit 

A kinematic fit has been carried out to improve the in- 
variant-mass and angular resolutions. This fit shifts the 
measured dK + missing mass (Fig. lb) to the nominal 
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Fig. 2. Upper: Dalitz plot of the events from the reaction 
pp^dK+K° at Q = 47.4 (left) and 104.7 MeV (right). The 
data are not background subtracted and not corrected for 
the ANKE acceptance and the detection efficiencies, and are 
binned with cell size 21 MeV 2 /c 4 x57MeV 2 /c 4 The lines denote 
the kinematically allowed region. Lower: Simulated Dafitz- 
plots inside the ANKE acceptance for phase-space distributed 
events (i.e. configuration [(KK) s d] s , cf. Sect. 2.5). 
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2.5 Acceptance correction 

In comparison to the data at T p = 2.65 GeV [21], the ex- 
cess energy for the higher beam energy is approximitely 
twice as large. As a consequence, the method of model- 
independent acceptance correction (using a five-dimen- 
sional acceptance matrix) can no longer be used, since 
the number of zero elements in the acceptance matrix be- 
comes too large. An alternative method has been devel- 
oped which is described as follows. 

In the close-to-threshold regime only a limited number 
of final states contribute. For the data analysis we have 
restricted ourselves to the lowest allowed partial waves, 
i.e. s-wave in the KK system accompanied by a p-wave 
of the deuteron with respect to the meson pair (aj (980)- 
channel) and p-wave KK production with an s-wave deu- 
teron (non-resonant channel). In the following we denote 
these two configurations by [(KK) s d] p and [{KK) p d] s . It 
has been shown that the lower energy data can be de- 
scribed by this ansatz for the dK + K° final state [21], 
where the square of the spin-averaged transition matrix 
element can be written as: 



\M\ 2 =C«q 2 + C*k 2 + C 1 (p-kf 
+C 2 (P ■ qf + C 3 (k ■ q) + C 4 (p • k)(p ■ q) 



(1) 



Here k is the deuteron momentum in the overall CMS, 
q denotes the K + momentum in the KK system, and p is 
the unit vector of the beam momentum. Only KK p- waves 
contribute to Cq and C 2 , only KK s- waves to Cq and C\ , 
and only s-p interference terms to C3 and C4. The coeffi- 
cients C'i can be determined from the data by fitting Eq.(l) 
to the measured dcr/dm K ^ and da/drridK as well as to 
the angular distributions d<r/d[cos (pk)], da/d[cos (pq)], 
dtr/d[cos (kq)] and dcr/d[cos (pt)] [19] (t represents the 
K + momentum measured in the overall CMS). It should 
be noted that a fit to the two-dimensional Dalitz plot does 
not provide additional information about the transition 
matrix, but would only yield three linear combinations of 
two of the coefficients C [19]. 

\M\ 2 gives the production probability of an event with 
certain kinematic parameters k and q relative to p. The 
corresponding differential acceptance of the spectrometer 
a(k,q,p) does not depend on the values of C, and can 
be determined using a large sample of simulated events, 
covering full phase space, which are tracked through a 
GEANT model of the setup [27]. Using the coefficients 
from Rcf. [21] as starting parameters, the simulations were 
carried out for different sets of the C leading to differ- 
ential distributions convoluted with the acceptance. For 
each choice of the C , the x 2 values have been calculated 
for the difference between simulated and measured distri- 
butions. Subsequently, the coefficients which describe the 
experimental data best have been determined by minimiz- 
ing x 2 with the MINUIT package [28]. The best fit result 
of this procedure is shown in Fig. 3 for two invariant-mass 
and four angular distributions (cf. Table 1 in Sect. 3.1 for 
numerical values). 
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Fig. 3. Best fit to non-acceptance-corrected data at T p = 2.83 
GeV. 



3 Cross sections of the reaction pp — > dK+K° 
at Q = 47.4 and 104.7 MeV 

3.1 Differential spectra 

With the best fit coefficients Ci one can simulate corre- 
sponding differential distributions at the target, track the 
events through the setup, and thus determine one-dimen- 
sional differential acceptances for e.g. the two invariant 
masses and four angles of Fig. 3. Using these acceptances, 
differential cross sections can be extracted from the data, 
and these are shown in Fig. 4. 

In order to verify the validity of the acceptance cor- 
rection method using the coefficients Ci, the same pro- 
cedure has been applied to the lower energy data. The 
results are shown in Fig. 5 as solid dots and are com- 
pared to our results published previously [21], analyzed 
using the acceptance matrix method (open circles). For a 
pp initial state all distributions must be forward-backward 
symmetric relative to the beam momentum; this feature 
has been exploited in Ref. [21], where differential cross 
sections as functions of | cos (pq)\ and | cos (pk)\ are pre- 
sented. These are shown in the lower left spectra of Fig. 5 
together with the mirrored distributions from the coeffi- 
cient method (squares), each scaled by 0.25 for better dis- 
tinction. In all cases good agreement between the model- 
independent matrix method [21] and the ansatz discussed 
here is observed. Note that the matrix method did not al- 
low us to extract the cos (pi) distribution (upper right in 
the Figure) from the 2.65 GeV data which is now possible 
with the coefficient ansatz. 

The best-fit coefficients Ci are presented in Table 1 
for both beam energies. All coefficients are given in units 
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Fig. 4. Angular and invariant mass distributions for T v — 2.83 
GeV. The dashed (dotted) line corresponds to KK production 
in a relative s- (p-)wave, the dash-dotted to the interference 
term, and the solid line is the sum of these contributions. The 
error bars show the statistical uncertainties only. The system- 
atic uncertainty for each bin is smaller than 10%; the over- 
all uncertainty from the luminosity determination is given in 
Sect. 2.1. 



of Cq. This is due to the fact that |A1| 2 from Eq.(l) is 
proportional to the differential cross sections, thus leaving 
one parameter undetermined in the fit. The errors of the 
d are obtained by varying each coefficient (allowing the 
others to change) such that the total x 2 increases by one. 

The parameters d from Eq.(l) can be directly related 
to the different partial waves [21]. Their contributions to 
the various observables are shown in_Fig. 4. The occur- 
rence of interference terms in the K°d invariant mass 
distribution is due to the choice of the kinematic vari- 
ables, i.e. relative momentum of the kaons and that of 
the deuteron with respect to the kaon pair ; Consequently, 
there is no interference term in the K + K° mass distri- 
bution. To get a distribution for the other invariant mass 
that is free of interferences one needs to switch to the 
K°d relative momentum and the K + momentum. Then, 
however, there will be an interference term in the K + K° 
mass distribution. The method how to construct the vari- 
able transformation is described in detail in Ref. [29]. 

Our fit reveals a strong dominance of the KK s-wave 
production rate (i.e. via the <2q (980) channel) for both 
beam energies: (95 ±4)% and (89 ±4)% at T p = 2.65 GeV 
and 2.83 GeV, respectively. 

The quality of the fit clearly supports the ansatz to 
include only the lowest partial waves in the data analysis. 
It should be noted that the growth of the amplitudes due 
to the centrifugal-barrier factor is taken care off by Eq.(l). 
An essential question is to understand the variation in the 



Fig. 5. Same as Fig. 4 for T v — 2.65 GeV but omitting the 
fitted partial-wave contributions. Full symbols denote the dif- 
ferential cross sections obtained by the method described in 
this paper; open symbols are the previously published model- 
independent results [21] where the error bars include statistical 
and systematic uncertainties. The cos (pt) spectrum has not 
been presented in Ref. [21]. See text for further details. 



parameters C3 and C4, see Table 1. As outlined above, 
these parameters emerge solely from an interference of the 
\(KK) s d] p with the [(KK) p d] s partial waves. Therefore, 
if there were a significant phase motion in one of these 
groups (e.g. due to the strong final state interaction in the 
ao channel), a variation with energy especially of C3 and 
C4 is expected. This point clearly calls for more theoretical 
investigations. 



3.2 Total cross sections 

Knowing the coefficients Cj, and thus the initial differen- 
tial distributions (Figs. 4 and 5), the total acceptance and 
the total cross sections can be evaluated. For the higher 
energy, a value of u(pp — ► dK + K°) = (190 ± 4 sta t ± 
39 sys t)nb is obtained. At the lower energy the extracted 
total cross section is in agreement with the previously pub- 
lished value of (38±2 st at±14 sys t) nb [21]. In both cases the 
errors include the statistical and systematic uncertainty 
from the luminosity determination. 

Figure 6 shows the measured total cross sections in 
comparison with the expected Q dependence of the cross 
section calculated with the transition matrix element of 
Eq.(l). After an angular integration the total pp — > dK + K° 
cross section is given by 



N 



a = 



2 6 7T 3 



k q 



4sTO 2 J V s S KK 

p Ami. 



\M\ 2 Askk , (2) 
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Table 1. Quality and results of the fit using Eq.(f). For the definition of N see Eq. 2. 
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Fig. 6. Total cross section of thepp— >dK + K° reaction as func- 
tion of the excess energy Q. The solid line is the result from 
Eq.(2) with the squared transition amplitude given by Eq.(3) 
and with iV=18 /A MeV -2 . The dashed line shows the energy 
dependence for three-body phase space. Note that the latter is 
forbidden by selection rules. 



where s and skk are the squared invariant energies of the 
initial pp and final K + K° systems, respectively. Here k 
and q are denned as before and are given explicitly as 

k 2 _ (f ~ s kk ~ m 2 d ) 2 - is KK m 2 d 



4s 



2 s K K - im K 
1 = a -■■ 



where rrid and ttik are the deuteron and kaon masses, i.e. 
we neglect the K + and K° mass difference. The angular- 
integrated squared transition amplitude \Ai\ 2 is given as 



\M\ 



(3) 



The normalization factor N has been determined for 
both energies and is quoted in Table 1. The errors of N 
include the systematic uncertainties of the total cross sec- 
tions. Note that the lines in Figs. 4 and 5 have been prop- 
erly scaled to the individual total cross sections. 

For illustration we show by the dashed line in Fig. 6 the 
result of Eq.(2) with a constant matrix element. This is a 
classical example when data can be well reproduced by a 
simple phase-space consideration although such a descrip- 
tion is invalid for this particular reaction since selection 
rules do not allow for a pure s-wave. 



4 Summary and conclusions 

Using the ANKE spectrometer at an internal target po- 
sition of COSY-Jiilich, we have searched for scalar KK 
production in the reaction pp — ► dK + K° at two excess 
energies Q = 47.4 and 104.7 MeV. Due to the excellent 
K + identification at ANKE, the detected events with co- 
incident K + d pairs exhibit little background. This can 
be subtracted using events outside the K° missing-mass 
peak^ After a kinematic fit to the K° mass, the invariant 
K + K° mass distribution has been obtained with an un- 
precedented resolution of better than 3 (10) MeV/c 2 for 
47.4 (104.7) MeV. 

Mass and angular distributions have been extracted 
from the data using an ansatz for the transition matrix 
element that includes the lowest allowed partial waves, 
i.e. an s-wave in the KK system accompanied by a p- 
wave of the deuteron with respect to the meson pair and 
p-wave KK production with an s-wave deuteron. All six 
coefficients that enter the spin-averaged matrix element 
have been obtained by a fit to the differential spectra. This 
fit reveals the dominance of KK production in a relative 
s-wave, (95 ± 4)% and (89 ± 4)% at 47.4 and 104.7 MeV, 
i.e. dominance of kaon-pair production via the a^(980) 
channel. 

The reaction pp — » dK + K° has been subject of several 
theoretical papers. For example, the authors of Ref. [30] 
point out that this reaction (and also pp — > dTr + r)) is ex- 
pected to be an additional source of information about 
the scalar sector. They account for the interactions of the 
mesons by using chiral unitary techniques, which dynam- 
ically generate the aj(980) resonance. In Ref. [31] total 
cross sections and differential spectra are calculated using 
a model in which the reaction pp^dK + K° is dominated 
by intermediate a^"(980) production. 

As discussed in Sect. 3.1, there is an energy dependence 
of the parameters C3 and C4 in the transition matrix el- 
ement (Eq.(l)) that is not yet understood and needs fur- 
ther theoretical study since it might indicate a final-state 
interaction in the ag channel. 
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